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Consider  the  problem  of  estimating  phase  and  decoding  data  symbols 
from  baseband  data.  Assume  the  phase  sequence  is  a random  walk  on 
the  circle  and  the  symbols  are  drawn  independently  from  an  equi- 
probable  alphabet  for  transmission  over  a perfectly  equalized  channel . 
A dynamic  programming  algorithm  (Viterbi  algorithm)  is  derived  for 
decoding  a maximum  a posteriori  (MAP)  phase-symbol  sequence  on  a 
finite  dimensional  phase-symbol  trellis.  y Simulation  results  for 
binary  and  8-ary  phase  modulated  (PM)  symbol  sets  are  presented. 


I.  Introduction 

Current  techniques  for  decoding  data  symbols 
transmitted  over  random  phase  channels  include  the 
decision-directed  phase-lock  loop  (DDPLL)  of  (1)  and 
the  decision-directed  stochastic  approximation  proced- 
ure of  [2).  The  potential  of  maximum  a posteriori 
(MAP)  sequence  estimation  for  improving  on  the  DDPLL 
has  been  recognized  in  13),  where  the  author  derives 
a path  metric  and  indicates  its  role  in  a forward 
dynamic  programming  algorithm  for  obtaining  MAP  phase- 
symbol  sequences.  However,  because  of  the  way  random 
phase  is  modelled  in  [3],  two  simplifying  assumptions 
must  be  made  in  order  to  implement  a tractable,  finite 
dimensional  algorithm. 


^(x)  = T. 


m=l  «=-<*>  (2r<j~)' 


—y~T  OX pi-- 


>:-:2’-(rrr-0  • 


We  may  think  of  {$,  } as  a modulo-2-  version  of  usual 
rai.dom  walk.  The  function  fcM(x)  is  useful  for  our 
discussion  of  PM  in  Section  TV. 

The  noise  sequence  is  assumed  to  be  a sequence  of 
complex  normal  random  variables: 


In  this  paper  we  observe  that  only  phase-modulo 
2tt  is  of  interest  in  data  communication  appl  lent  Ions . 
This  motivates  us  to  wrap  the  phase  around  the  circle, 
so  to  speak,  and  obtain  a folded  normal  model  for 
transition  probabilities  in  random  walk  on  the  circle. 
It  is  then  straight-forward  to  pose  a MAP  phase-symbol 
sequence  estimation  problem  as  in  [5]  and  [6).  The 
basic  idea  is  to  discretize  the  phase  space  [-71,71)  to 
a finite  dimensional  grid  and  to  use  a dynamic  pro- 
gramming algorithm  (Viterbi  algorithm)  to  keep  track 
of  survivor  phase-symbol  sequences  that  can  ultimately 
approximate  the  MAP  sequence. 


II.  Model  for  Symbol  Trnnsmi ss ion  over  a Random  Walk 
Phase  Channel 

As  our  model  for  symbol  transmission  over  a per- 
fectly equalized  channel  we  assume 


"k  = "k  + 1vk 


IL  v(  V(k..> 


zk  = ‘k  e 


+ n.  , k = 1 , 2 , . . . 
k 


Here  (z^)  *s  a sequence  of  baseband  measurements,  fi^) 
is  an  Independent  sequence  of  M-.iry  complex  symbols 
drawn  from  an  oqul probab  1 e alphabet,  {.f»  I Is  .1  random 
walk  phase  sequence,  and  f n^ I Is  an  additive  noise 
sequence.  The  sequences  (I  ^ I , I I , and  1 11  I are  (n 
dependent  of  each  other. 

Tin*  phase  is  modelled  as  a random  walk  on  the 
circle  taking  values  in  [-77, it).  This  sequence  is 
Markov  with  transition  probabilities  characterized  by 
the  following  folded-normal  law  [4],  [5]: 
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Here  u :N(0,o  ) indicates  u,  is  normal  with  mean  zero 
k n_  R 

and  variance  0“.  It  follows  that  the  joint  density  of 
^ K K 

the  observations  {z^j*  phases  { ; ^ ^ , and  symbols 
{1^}^  may  be  written 

V?-  (V?>  ■ 


M"  ^Vv  ’V*i<vwf<v 

where  f (4>  ) is  the  marginal  density  of  $ and  where 

2 ^ ^2-1 
N_  (x,o  ) denotes  the  complex  normal  density  (2na  ) 

7V  n 2 2 n 

expi-[z^-x|  /2oM.  The  Joint  density  *1 1 (4)  is  actu- 
al I V a mixed  dens  1 1 v-pi'ob.ibl  lily  miss  I unci  inn  defined 
on  K^x|  n,n)^  >.  1^  where  K Is  I he  real  line  and  I Is 
the  M -d i mens  I on.  1 1 discrete-set  of  source  symbols. 


ITT.  A Viterbi  Algorithm  for  MAP  Phase  and  Symbol 
Sequence  Decoding 

Consider  the  following  maximization  problem  with 
respect  to  the  phase  sequence  and  the  symbol 


sequence 


!“x  „ f<lzk)l’  (Vo  ‘Vl*  • <5) 

Maximization  of  this  joint  density  function  is  equiva- 
lent to  maximization  of  the  a posteriori  density 

fd^k^l*  80  we  call  the  maximizing 

sequences  {1^}J  the  MAP  phase-symbol  sequences. 

If  instead  we  maximize  the  in  f(.,  . ,.)  and  ignore  un- 
interesting constants  we  have  the  problem 


u .K,  (1  >K 

l'k  r 1 k i 

inhere  is  defined  recursively  in  terms  of  the  path 


K J*k  2 -1 

fK  . * Nzk  . V « S1(R(VVl))  (12) 

where  K(x)  Is  the  difference  between  x and  IxJ.  Thus 

K 

we  may  maximize  (12)  with  respect  to  {|J>  } , and  then 

* K k 1 

find  {0  1.  from  (11).  In  a sense,  we  are  decoding 

K * * 

• as  if  it  had  a "density"  given  by  M gj(K(*)) 
and  simply  using  the  algorithm  of  [ 5 J and  (6).  Wien 
the  overlap  between  2ti/M  translates  of  (^“^  j)  is 
negligible,  th ;n  KM(*k-'f,k_l ) 4 M 1 Ki (R (!^k_'f'k-l > > and 
we  may  write  the  problem  of  maximizing  £n  f^  as  in  (6) 
and  (7)  with  replaced  by  p': 


+ eM(Vl'k-l) 
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In  this  form  the  maximization  problem  may  be  efficient- 
ly solved  by  using  a dynamic  programming  algorithm  (or 
Viterbi  Algorithm)  on  a finite-dimensional  grid  of 
phase-symbol  pairs.  The  reader  is  referred  to  [5]  and 
[6]  for  details. 

IV An  Optimum  Decoding  Algorithm  for  Phase-Coded 
Symbols 

. c t r j2n (m-l)/M,M 

Suppose  l = (e  'm=l  corrcsPoncling  to 

equally-spaced  PM  symbols  with  unit  energy.  Write 

iK 


1,2,... 


*k  + 9k 


and  consider  the  following  MAP  sequence  estimation 
problem: 

“x  K f({zk}?>  {Vr  {V^  (9) 

(Vi*  <vf 

Assuming  {t^l  to  be  defined  on  ( — tt , tt ) one  can  show 
that  the  joint  density  in  (9)  may  be  written 

K 2 -1 

fK  = k!1  \ <e  ’°n)  M «l(VVf(Vek-l»' 

(10) 

* K “ K 

Call  { i|'k ) j and  {9^}^  the  MAP  sequences  that  jointly 
maximize  f^.  These  sequences  enter  jointly  only  in 
the  gj(*)  term  on  the  right  hand  side  of  (10).  It 
follows  then  that  f^  is  minimized  by  choosing. 

\ - Vi  + 'VVi1  <1J> 

win* re  |*1  denotes  the  closest  value  of  (m-I)  2it/M  to 

i|»k  -i|»k  j and  0^  is  Interpreted  modiiln-?n . Ii  follows 

f,  1n.1v  he  rewritten  as 
k 7 


The  function  g^  is  the  convolution  of  the  folded  phase 
density  with  the  impulsive  density  of  the  data. 

V.  Results 

The  phase  space  [-tt,tt)  has  been  discretized  to 
48  phase  values  and  a Viterbi  algorithm  has  been  pro- 
grammed to  solve  (9)  along  the  lines  outlined  in  (11)- 
(13).  Initial  phase  acquisition  has  been  achieved  by 
sending  a 100-symbol  preamble  followed  by  900  data 
symbols.  This  procedure  has  been  repeated  from  10  to 
20  times  to  obtain  the  performance  results  of  Figs.  1 
and  2.  When  counting  symbol  errors,  bursts  have  been 
ignored.  Shown  in  Fig.  1 are  binary  symbolling  re- 
sults for  0^  = 0.01  rad^  (o  =5.7°)  and  SNR  = 10  log,_ 

2 w w 10 

l/2o  ranging  from  4 to  lOdB.  The  results  for  binary 
n 

orthogonal  symboling  are  of  no  inherent  interest  in 
their  own  right.  They  are  presented  simply  to  vali- 
date the  simulation.  The  results  for  binary  orthogo- 
nal symboling  are  interesting  because  they  indicate 
performance  nearly  equivalent  to  a receiver  that  has 
perfect  phase  coherence.  The  circles  of  Fig.  2 repre- 
sent error  probabilities  for  8-ARY  PM  symboling  when 

2 2 ^ 

SNR  ranges  from  16-19dB  and  (a  0 )^  remains  fixed  at 
-3  2 n w 

4.4  x 10  rad  . The  triangles  represent  error  proba- 
bilities for  the  markedly  simpler  decision-directed 

algorithm  of  [21.  For  pure  PM  and  moderate  values  of 

2 

the  phase  variance  parameter  o , there  seems  to  be 

w 

little  to  recommend  the  algorithm  of  Sections  III  and 
IV  over  the  algorithm  of  [2).  A similar  conclusion 
with  regard  to  the  DDPI.L  was  reached  in  [31.  Tills 
conclusion  is  changed  for  AM-PM  symbol  constol lat ion 
such  as  lb-QASK,  4A-44*,  etc. 
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